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Using a formulation of the post-Newtonian expansion in terms of Feynman graphs, we discuss 
how various tests of General Relativity (GR) can be translated into measurement of the three- 
and four-graviton vertices. In problems involving only the conservative dynamics of a system, a 
deviation of the three-graviton vertex from the GR prediction is equivalent, to lowest order, to 
the introduction of the parameter /3ppn in the parametrized post-Newtonian formalism, and its 
strongest bound comes from lunar laser ranging, which measures it at the 0.02% level. Deviation of 
the three-graviton vertex from the GR prediction, however, also affects the radiative sector of the 
theory. We show that the timing of the Hulse- Taylor binary pulsar provides a bound on the deviation 
of the three-graviton vertex from the GR prediction at the 0.1% level. For coalescing binaries at 
interferometers we find that, because of degeneracies with other parameters in the template such 
as mass and spin, the effects of modified three- and four-graviton vertices is just to induce an error 
in the determination of these parameters and, at least in the restricted PN approximation, it is not 
possible to use coalescing binaries for constraining deviations of the vertices from the GR prediction. 
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I. INTRODUCTION 

Binary pulsars, such as the Hulse- Taylor [1] and the 
double pulsar [2], are wonderful laboratories for testing 
General Relativity (GR). They have given the first ex- 
perimental confirmation of the existence of gravitational 
radiation [3, 4], provide stringent tests of GR and allow 
for comparison with alternative theories of gravity, such 
as scalar-tensor theories [5-14] (see [15-17] for reviews). 
Another very sensitive probe of the non-linearities of GR 
is given by the gravitational wave (GW) emission during 
the last stages of the coalescence of compact binary sys- 
tems made of black holes and/or neutron stars, which is 
one of the most promising signals for GW interferome- 
ters such as LIGO and Virgo, especially in their advanced 
stage, and for the space interferometer LISA. Various in- 
vestigations have been devoted to the possibility of using 
the observation of coalescing binaries at GW interferom- 
eters to probe non-linear aspects of GR [9, 10, 18-24]. 

Compact binary systems probe both the radiative sec- 
tor of the theory, through the emission of gravitational 
radiation, and the non-linearities intrinsic to GR which 
are already present in the conservative part of the La- 
grangian. In a field-theoretical language, these non- 
linearities can be traced to the non-abelian vertices of 
the theory, such as the three- and four-graviton vertices. 
It is therefore natural to ask whether from binary pul- 
sars or from future observations of coalescing binaries at 
interferometers one can extract a measurement of these 
vertices, much in the same spirit in which the triple and 
quartic gauge boson couplings have been measured at 
LEP2 and at the Tevatron [25-28]. 

In this paper we tackle this question. The organiza- 
tion of the paper is as follows. In Section II we dis- 
cuss how to "tag" the contribution of the three- and 



four-graviton vertices to various observables in a con- 
sistent and gauge-invariant manner, and we compare it 
with other approaches, such as the parametrized post- 
Newtonian (PPN) formalism [16]. In particular, we find 
that the introduction of a modified three-graviton vertex 
corresponds - in the conservative sector of the theory 
and at first Post-Newtonian order (1PN) - to the intro- 
duction of a value for the PPN parameter /3ppn different 
from the value /?ppn = 1 of GR. However, a modified 
three-graviton vertex also affects the radiative sector of 
the theory, which is not the case for the PPN parameter 
/3ppn- We also discuss subtle issues related to the possi- 
ble breaking of gauge invariance which takes place when 
one modifies the vertices of the theory. In Section III we 
present our computations with modified vertices, and in 
Section IV we compare these computations with exper- 
imental results obtained from the timing of binary pul- 
sars and with what can be expected from the detection of 
gravitational waves (GWs) at ground-based interferom- 
eters or with the space interferometer LISA. Section V 
contains our conclusions. 



II. TAGGING THE THREE- AND 
FOUR-GRAVITON VERTICES 

Our aim is to quantify how well the non-linearities of 
GR can be tested by various existing or planned experi- 
ments/observations. Historically, there have been several 
approaches to this problem and, basically, one can iden- 
tify two complementary strategies. The first is to develop 
a purely phenomenological approach in which deviations 
from GR are expressed in terms of a number of parame- 
ters, without inquiring at first whether such a deforma- 
tion of GR can emerge from a fundamental theory. An 
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example of such an approach is the parametrized post- 
Newtonian (PPN) formalism. In its simpler version, it 
consists of writing the 1PN metric generated by a source, 
treated as a perfect fluid with density p(x) and velocity 
field i>(x), in the form 
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and the standard PPN gauge has been used [16, 29] 
(we use units c = 1). General Relativity corresponds 
to /?ppn = 1 and 7ppn = 1- More phenomenological 
parameters can be introduced by working at higher PN 
orders, see [16]. One then investigates how deviations of 
/?ppn and 7ppn from their GR values affect various ex- 
periments. Writing /?ppn = 1 + (3 and 7ppn = 1 + 7, the 
best current limits (at 68% c.l.) are 



7= (2.1 + 2.3) x 10" 



(6) 



from the Doppler tracking of the Cassini spacecraft, and 

4/3-7 = (4.4 + 4.5) x 10" 4 (?) 

from lunar laser ranging. This bound comes from the 
Nordtvedt effect, i.e. from the fact that, in a theory with 
(3 and 7 generic, the weak equivalence principle is vio- 
lated and the Earth and the Moon can fall toward the 
Sun with different accelerations, which depend on their 
gravitational self-energy. The effect is studied by moni- 
toring the Earth-Moon distance with lunar laser ranging. 
The perihelion shift of Mercury gives instead the bound 
< 3 x 1(T 3 [16, 30]. 
To get these bounds, we do not need to know the funda- 
mental theory that gives rise to values of /3ppn and 7ppn 
that differ from their GR values. However, it is of course 
interesting to see that consistent field theories exist that 
give rise to values of /3ppn and 7ppn different from one. 
For instance, a Brans-Dicke theory with parameter wbd 
gives /3ppn = 1 and7ppN = (1+^ B d)/(2+uibd), with the 
GR value 7ppn = 1 recovered for wbd — > 00, while more 
general tensor-scalar theories can produce both 7ppn 7^ 1 
and /3ppn 7^ 1- However, in the PPN approach, one can 
also explore other possibilities, such as PPN parameters 
that correspond to preferred-frame effects or to viola- 
tion of the conservation of total momentum, which are 
not necessarily well-motivated in terms of current field- 
theoretical ideas on possible extensions or UV comple- 
tions of GR. It is also important to observe that the pa- 
rameters /3ppn and 7ppn are gauge-invariant, and there- 
fore observables, because they have been defined with 



respect to a specific gauge, namely the standard PPN 
gauge in which the metric takes the form (l)-(3). 

A second, complementary, approach to the problem 
is to study a specific class of field-theoretical extensions 
of GR. A typical well-motivated example is provided by 
multiscalar-tensor theories. These have been studied in 
detail and compared with experimental tests of relativis- 
tic gravity in rcfs. [6, 9-14]. This approach has the ad- 
vantage that one is testing a specific and well-defined 
fundamental theory. On the other hand, an experimen- 
tal bound on the parameters of a given scalar-tensor the- 
ory, such as for instance the bound wbd > 40000 on the 
parameter wbd of Brans-Dicke theory obtained from the 
tracking of the Cassini spacecraft [31] is, strictly speak- 
ing, only a statement about that particular extension of 
GR and not about GR itself. 

In this paper we quantify how well GR performs with 
respect to experiments of relativistic gravity by studying 
how much these experiments constrain the values of the 
non-abelian vertices of the theory, in particular the three- 
graviton vertex and the four-graviton vertex. We proceed 
as follows. After choosing a gauge (the De Donder gauge, 
corresponding to harmonic coordinates) we multiply the 
three-graviton vertex by a factor (1 + (3 3 ) and the four- 
graviton vertex by a factor (1 + ^4), with constants fa 
and Pi. For fa — $4 — we recover GR. Observe that, 
since fa and fa are defined with respect to a given gauge 
choice, they are gauge-invariant by definition. This is in 
fact the same logic used to define in a gauge-invariant 
manner the PPN parameters /3ppn and 7ppn- 

We then use a Feynman diagram approach to com- 
pute the modifications induced by fa and fa on various 
observables in classical GR. Diagrammatic approaches 
and field-theoretical methods have been in use in clas- 
sical GR for a long time, see e.g. [13, 32, 33]. We make 
use of the effective field theory formulation proposed 
in [34] , which provides a clean and systematic separation 
of the effects that depend on (model-dependent) short- 
distance physics from long-wavelength gravitational dy- 
namics, and in the non-relativistic limit (after perform- 
ing a multipole expansion) has manifest power counting 
in the typical velocity v of the source. 

In the next section we see how these deformed vertices 
give additional terms in the PN effective Lagrangian. 
In particular we find that, in the conservative sector of 
the theory, the introduction of fa is phcnomcnologically 
equivalent, at 1PN level, to the introduction of a non- 
trivial value of /3ppN given by /3ppn = 1 + fa- However, 
fa also affects the radiative sector of the theory, i.e. the 
Lagrangian describing the interaction between the mat- 
ter fields and the gravitons radiated at infinity. 

Before entering into the technical aspects, however, let 
us clarify the meaning of the introduction of fa and (3±. 
In ordinary GR, with fa^ = 0, coordinate transforma- 
tion invariance ensures that the negative norm states de- 
couple. After gauge fixing (in the De Donder gauge for 
instance), the kinetic terms for all of the ten components 
of the metric are invertible, but four of them have the 
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wrong sign, i.e. they give rise to negative norm states. 
In the De Donder gauge the six positive-norm states that 
diagonalize the kinetic term are 

hij = hij + -Sij(h o — S lm hi m ) , (8) 

while the four "wrong-sign" components are given by the 
spatial vector h oi and by the scalar 

h N = h 00 - S lm h lm . (9) 

In standard GR the existence of these negative-norm 
states do not create difficulties because they are cou- 
pled to four integrals of motion (energy, and the three 
components of angular momentum), so they cannot be 
produced. In contrast, the remaining six "healthy" com- 
ponents couple to the source multipole moments. Af- 
ter complete gauge fixing one finds that among the six 
positive-norm states, four obey Poisson-like equations, so 
they do not radiate (even though they are non-radiative 
physical degrees of freedom), while the remaining two are 
the radiative degrees of freedom representing GW's [35]. 

Allowing ft^O has the effect that the negative norm 
state Jin now couples, already at lowest order, to a non 
conserved quantity, namely to a combination of the New- 
tonian kinetic and potential energy of the binary system. 
This means that, in general, a modification of GR in 
which we just change the strength of the three-graviton 
vertex cannot be taken as a fundamental field theory, nei- 
ther at the quantum level, nor even at the classical level, 
since the negative-norm state contributes to the classical 
radiated power (a related concern is that, for 3 ^ 0, the 
energy-momentum tensor is in general not conserved). 
A consistent classical and quantum field theory could in 
principle emerge from a simultaneous modification of all 
the vertices of the theory, such as the three-, four- and 
higher-order graviton vertices, together with a related 
modification of the graviton-matter couplings. As a triv- 
ial example, an overall rescaling of the gauge coupling in a 
Yang-Mills theory, or of Newton's constant in GR, results 
in a combined modification of all the vertices, but obvi- 
ously introduces no pathology. Anyway, our approach 
to the problem is purely phenomenological. We intro- 
duce /?3 and 04 simply as "tags" that allow us to track 
the contribution of the three- and four-graviton vertices 
throughout the computations. As long as | /?3 1 <C 1 and 
|/34 1 <C 1, the corrections that they induce to the radiated 
power are small compared to the standard GR result, so 
the total radiated power is given by the GR result plus 
a small correction, and in particular the total radiation 
emitted is positive. At this phenomenological level the 
introduction of modified vertices is therefore acceptable, 
and provides a simple and, most importantly, gauge in- 
variant manner of quantifying how well different obser- 
vations constrain the non-linear sector of GR, in a way 
which is intrinsic to GR itself, without reference to any 
other specific field theory. 

In this sense, our approach is close in spirit to the 
phenomenological PPN approach, and can be seen as an 



extension of it where the radiative sector of the theory 
is also modified. Another approach which is related to 
ours is the one proposed in ref. [23, 24]. They consider 
the phase of the GW emitted during the coalescence of 
compact binaries, which up to 3.5PN has the form 

7 

*(/) = 27r/i c - $ c + Y}4>k + ln/]/ (fc - 5)/3 , (10) 

fe=0 

where / is the GW frequency, and t c and $ c are the time 
and the phase at merger. The seven non-zero coefficients 
tpk with k = 0, 2, 3, . . . 7 and the two non-zero coefficients 
tpki with k — 5,6 are known from the PN expansion, in 
terms of the two masses mi and mi. In ref. [23, 24] they 
study how the template is affected if these coefficients 
are allowed to vary, so that two of them, the OPN coeffi- 
cient tpo and the 1PN coefficient ip2, are used to fix the 
masses toi and m 2 of the two stars, while varying any 
of the remaining coefficients with k > 3 provides a test 
of GR. In the case of coalescing binaries at interferome- 
ters, our introduction of 03 and 04 is a particular case of 
a more general analysis in which one treats the quanti- 
ties ipk and ipki as free parameters, but it has a sharper 
field-theoretical meaning since 03 and 04 measure the de- 
viation of the three- and four-graviton vertices from the 
GR prediction. For the same reason, we are also able to 
compare the effect of 03 on the waveform of coalescing 
binaries with its effect on binary pulsar timing and on 
solar system experiments, while in the phenomenological 
approach in which the parameters ^fc and ipki of the GW 
phase are taken as free parameters, a modification of the 
waveform of coalescing binaries cannot be related to a 
modification of the binary pulsar timing formula. 

Another issue is whether a modification of the ver- 
tices of this form (typically with 03 , 04, etc. not inde- 
pendent, but related to each other by some consistency 
conditions) could emerge from a plausible and consistent 
extension of GR. Actually, a typical UV completion of 
GR at an energy scale A will rather generate corrections 
to the vertices that are suppressed by inverse powers of 
A, so it would give rise to an energy dependent 03, e.g. 

03 = E 2 /A 2 , which furthermore, at the energy scales 
that we are considering and for any sensible choice of 
A, would be utterly negligible. Still, let us remark that 
this kind of behavior is not a theorem. It assumes the 
UV-IR decoupling typical of effective field theories, and 
one can exhibit counterexamples. For instance, in non- 
commutative Yang-Mills theories there is a UV-IR mix- 
ing, such that low-energy processes receive contributions 
from loops where very massive particles are running, and 
these contributions are independent of the mass of these 
particles [36] . Anyway, again our aim here is not to test 
any given consistent extension of GR, but rather pro- 
vide a simple and phenomcnologically consistent way of 
quantifying how well various experiments can test the 
non-linearities of GR, and quantify how the results of 
different experiments compare among themselves. 

It is also interesting to observe that, even when 03 and 

04 are non-zero, the graviton remains massless at the 
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classical level, since fa and fa affect interaction terms, 
but not the kinetic term. The breaking of diffcomor- 
phism invariance induced by fa and fa could in principle 
generate a graviton mass at the one-loop level. However, 
even if we are using the language of quantum field theory, 
in the end we are only interested in the classical theory, 
since quantum loops are suppressed by powers of fi/L, 
where L is the angular momentum of the system, so they 
are completely negligible for a macroscopic system. [55] 



III. EFFECTIVE LAGRANGIAN FROM A 
MODIFIED THREE-GRAVITON VERTEX 

To perform our computations we use the effective field 
theory formulation proposed in ref. [34]. Computations 
of the conservative dynamics at 2PN level have been per- 
formed using this effective field theory technique [37, 38] 
and the results are in agreement with the classic 2PN 
results of refs. [39, 40], while the full 3PN result for non- 
spinning particles to date has only been obtained with the 
standard PN formalism using dimensional regularization 
[32, 33] (see ref. [46], or Chapter 5 of ref. [17], for a peda- 
gogical introduction to the PN expansion and for a more 
complete list of references). Spin-spin contributions at 
3PN order in the conservative two-body dynamics have 
recently been computed both with the effective field the- 
ory techniques [41, 42] and with the ADM Hamiltonian 
formalism [43] (see also [44, 45] for other applications of 
the EFT technique related to gravitational radiation). 

In the formalism of ref. [34], after integrating out 
length-scales shorter than the size of the compact ob- 
jects, the action becomes 

S = Seh + S pp , (11) 

where Seh is the Einstein- Hilbcrt action and 



S n 



(12) 



is the point particle action. Here a = 1, 2 labels the two 
bodies in the binary system and dr a — y/ 'g^^x^dxadx^. 
One then observes that, in the binary problem, the gravi- 
tons appearing in a Feynman diagram can be divided into 
two classes: the forces between the two bodies with rela- 
tive distance r and relative speed v are mediated by gravi- 
tons whose momentum scales as (k° ~ v/r, |k| ~ 1/r). 
These are called "potential gravitons" and are off-shell, 
so they can only appear in internal lines. The gravitons 
radiated to infinity rather have (k° ~ v/r, |k| ~ v/r). 
One then writes =_r)^v + h^v and separates into 
two parts, hfj, v {x) — h^ v {x) + H tllJ {x) with h^ u {x) de- 
scribing the radiation gravitons and (x) the potential 
gravitons. One fixes the de Donder gauge and, expand- 
ing the action in powers of h^ v and Hy Siflu (xo), one can 
read off the propagators and the vertices, and write down 
the Feynman rules of the theory. Then, using standard 
methods from quantum field theory, one can construct an 
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FIG. 1: The diagrams that give the terms proportional to 
G 2 N in the 1PN Lagrangian. Dashed lines denote potential 
gravitons, solid lines the point-like sources. 



effective Lagrangian that, used at tree level, reproduces 
the amplitudes computed with the Feynman graphs. For 
a classical system, whose angular momentum L h, 
only tree graphs contribute, and reproduce the classical 
Lagrangian that is usually derived from GR using the 
PN expansion. For instance, the conservative dynamics 
of the two-body problem, at 1PN level, is given by the 
Einstein Tnfcld - Hoffmann Lagrangian, 
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In the language of the effective field theory of Ref. [34], 
this result is obtained from Feynman diagrams involving 
the exchange of potential gravitons. In particular, the 
terms linear in Gjv in eq. (13) are obtained from a sin- 
gle exchange of potential gravitons between two matter 
lines, see Fig. 4 of Ref. [34] , while the term proportional 
to G 2 N is obtained from the sum of the two graphs shown 
in Fig. 1. In the derivation of the conservative 1PN La- 
grangian, the three-graviton vertex only enters through 
the graph in Fig. la. Multiplying this vertex by a fac- 
tor (1 + fa) and repeating the same computation as in 
ref. [34] we get the additional contribution to the conser- 
vative part of the Lagrangian 
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G 2 N mim 2 (mi + m 2 ) 



(14) 



Comparing this result with the Lagrangian whose 
equations of motion are the same as the equations of 
motion of a test particle in the PPN metric (l)-(3) (see 
eq. (6.80) of ref. [29]) we find that, to 1PN order and as 
far as the conservative dynamics is concerned, the intro- 
duction of fa gives rise to a PPN theory with j3 — l+fa, 
i.e. j9 = fa, while 7 = 1 as in GR. Therefore the 
bound on (3 from the perihelion of Mercury translates 
into \fa\ < 3 • 10~ 3 , while eq. (7) translates into the 
bound (at 68% c.l.) 



\fa\ < 2-10- 



(15) 



This bound reflects the fact that fa, just as the PPN pa- 
rameter /?, violates the weak equivalence principle. The 



introduction of ft, however, also affects the radiative sec- 
tor of the theory, something that is not modeled in the 
phenomenological PPN framework since the latter by def- 
inition is only concerned with the motion of test masses 
in a deformed metric, and therefore only modifies the 
conservative part of the dynamics. Note however that, 
in the framework of multiscalar-tensor theories, the ex- 
tension of the PPN formalism introduced in Ref [11] al- 
lows for a consistent treatment of both the conservative 
dynamics (including the case of strongly self-gravitating 
bodies) and of radiative effects. 

It is clearly interesting to see what bounds on ft can 
be obtained from experiments that probe the radiative 
sector of GR, such as the timing of binary pulsars or the 
observation of the coalescence of compact binaries at in- 
terferometers. The effective Lagrangian describing the 
interaction of the binary system with radiation gravitons 
is obtained by computing the three graphs in Fig. 2 (cor- 
responding to Fig. 6 of ref. [34]), and the introduction of 
ft affects the HHh vertex in Fig. 2c. 

Computing these graphs as in ref. [34], but with our 
modified three-graviton vertex, we find 

•Crad = f,»r [QijRoiOj + qRoiOi + ft^V/loo + Z lJ hij)] , 
ZMpi 

where Qij is the quadrupole moment of the source and 
we define 



V(r) 
Z ij (r) 



(17) 

(18) 
(19) 



where r = xi — X2. The term QijRoiOj in eq. (16) is the 
usual quadrupole interaction. The second term, qRom , is 
non-radiating when ft — 0, but we will sec that for ft 
it contributes to the radiated power when the orbit 
is non-circular. The last two terms in eq. (16) are the 
explicit ft-dependent terms induced by the modification 
of the three-graviton vertex in Fig. 2c. 

Finally, we have omitted a term where /loo is coupled 
to the conserved energy (at this order) and hoi is coupled 
to the conserved angular momentum, since these terms 
do not generate gravitational radiation. 

The back-reaction of GW emission on the source can 
be computed as usual from the energy balance equation 
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where Pgw is the power radiated in GWs and E the or- 
bital energy of the system. To obtain the expression for 
the radiated power for ft ^ we cannot simply use the 
quadrupole formula of GR, since the introduction of ft 
generates new contributions. To take them into account 
we proceed as in ref. [34], and we compute the imagi- 
nary part of the graph shown in Fig. 3. The vertices of 
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FIG. 2: The diagrams that contribute to the matter-radiation 
Lagrangian. Dashed lines denote potential gravitons, wiggly 
lines radiation gravitons, and solid lines the point-like sources. 



the graph can be read from eq. (16). When ft = the 
only relevant vertex comes from QijRaiQj, and the com- 
putation of the imaginary part gives back the Einstein 
quadrupole formula [34]. In our case we have various 
possible vertices, and we must compute the imaginary 
part of 



£ /^^/^(^(^(^(tx)^^)), (21) 
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matter variables and — (RoiOj,RoiOj,Sijhoo,hij) on 



depends on the 

the gravitational field. When both vertices of the dia- 
gram in Fig. 3 are proportional to the quadrupole, one 
obtains the usual GR result 
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g (Qij Qij) 
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as already found in [34] . Computing the other contribu- 
tions we find that the terms Pq ? and P qq vanish identi- 
cally. In fact, the Qq and qQ graphs vanish because Qij 
is traceless, while the qq graph vanishes because SijSki 
gives zero when contracted SikSji + SuSjk — which 
is the tensor that comes out from the two-point function 
(RoiOj Rokoi) ■ The QV, qZ and VZ graphs vanish for 
similar reasons, so the only relevant contributions come 
from the QZ and qV graphs, and we find 



and 



Pqz — —ZfoGNiQijZij) 



P g v = -6f3 3 G N (qV) 



(23) 



(24) 



As for the VV and ZZ graphs, they give a contribution 
that, from the point of view of the multipole expansion, 
is of the same order as the quadrupole radiation but pro- 
portional to ft 2 , and can be neglected. 

We can now use these results to perform the compari- 
son with binary pulsars and with interferometers. 



IV. COMPARISON WITH EXPERIMENT 

As we already saw in eq. (15), solar system experi- 
ments, and in particular lunar laser ranging, give the 
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FIG. 3: The self-energy diagrams whose imaginary part gives 
the radiated power. The wiggly line can refer either to hij 
or to hoo, and the vertices to any of the four IijS^ with 
a — 1, .... 4, see the text. 



bound < 2-10~ 4 . In this section we study the bounds 
on /?3 that can be obtained from binary pulsars and from 
the detection of coalescing binaries at interferometers. 



The term proportional to /(e) is the standard GR re- 
sult [47], while the term proportional to g(e) is the extra 
contribution due to ^3. 

In order to measure fa from the dynamics of binary 
pulsars, we must also determine the dependence on (3$ 
of the periastron shift u and of the Einstein time delay 
7_e, since these two observables are used to determine 
the masses of the two compact stars. In particular the 
periastron shift fixes the total mass M of the system, 
while the Einstein time delay 7^ measures a different 
combination of masses, see e.g. eqs. (6.56) and (6.93) 
of ref. [17]. Using the conservative Lagrangian with the 
modification (14) and repeating the standard textbook 
computation of the periastron shift u, we find that the 
value LOp 3 computed in a theory with (3% ^ is related to 
the GR value wgr by 



A. Binary pulsars 



I 1 03 , 



(31) 



Since ^3 modifies the emitted power already at Newto- 
nian order, the energy of the orbit in eq. (20) can now be 
directly computed using the Keplerian equations of mo- 
tion. We see that this test of GR is conceptually different 
from the tests based on solar system experiments. The 
latter only probe the conservative part of the Lagrangian, 
i.e. the /33-dependent term given in eq. (14), while binary 
pulsars are sensitive to the /3 3 dependence given in the 
radiation Lagrangian (16) (even if the effect of f3s on the 
conservative dynamics will also enter, through the deter- 
mination of the masses of the stars from the periastron 
shift, see below). 

Using the Keplerian equations of motion for an elliptic 
orbit of eccentricity e we get 



Pqq 

Pqz 



32G%fi 2 M 3 
5a 5 (l - e 2 ) 7 /2 
32G%fi 2 M 3 



1 



73 

— e" H e 

24 96 



37 4 







P 



qV 



where M 



-0: 



5a 5 (l - e 2 ) 7 /2 



175 
~24 t 



85 
96* 



5a 5 (l-e 2 ) 7 / 2 ^16 



64 



(25) 
,(26) 
(27) 



mi + ttt-2 is the total mass, [i = m±m,2/M 
is the reduced mass, and we will also use the notation 
v = m\milM' 2 for the symmetric mass ratio. From the 
energy balance equation we then get the evolution of the 
orbital period P/,, 



n 

where 







-8/3 



[/(e) + 039(e)], (28) 



/(e) 



3(e) 



(1 - e 2 ) 7 / 2 



(1 - e 2 ) 7 / 2 



73 2 

24 6 

335 



37 
96 f 



155 
192 f 



(29) 
(30) 



while the Einstein time delay is unchanged because it 
is not affected by the post-Keplerian parameters. So, if 
/?3 7^ 0, the true value of the total mass M of the binary 
system, that enters in eq. (28), is not the one that would 
be inferred from the periastron shift using the predictions 
of GR, but rather we get 



! 1 + y ) Ugh 



(32) 



Similarly, using eq. (6.93) of ref. [17], for the symmetric 
mass ratio v we get 



v (3 3 = (1 + W03)VGK: 



where 



and 



k y/l + 4k - 2 



3 V1 + 4k (1+4k) 1 /2_(i + k ) 



.= 2(|) V3 (G ,M GR )-V3 



(33) 



(34) 



(35) 



Putting everything together and keeping only the linear 
order in fa we finally find that the ratio between the value 
of Pb computed at /3 3 ^ and the value of computed 
in GR is 



P 



(ft.) 



I + 039(e), 



where 



-/ \ 9(e) 5 
9{e) = W) + 6 +W - 



(36) 



(37) 



Observe that the term g(e)/ f(e) comes from the effect of 
/3 3 on the radiative sector of the theory, while the term 
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(5/6) + w comes from the effect of fa on the conser- 
vative sector, i.e. on the mass determination. Insert- 
ing the numerical values for the Hulse- Taylor pulsar, we 
get g(e) ~ 3.21. Note that g(e)/f(e) ~ 2.38, so 5(e) is 
dominated by the effect of fa on the radiative sector of 
the theory. For this binary pulsar, after correcting for 
the Doppler shift due to the relative velocity between 
us and the pulsar induced by the differential rotation of 
the Galaxy, the ratio between the observed value p° bs 
and the GR prediction P b GR is p°bs/j>GR = i.Q013(21). 
Interpreting this as a measurement of fa we finally get 
3.21/? 3 = 0.0013(21), i.e. 



fa = (4.0 ±6.4) • 10" 



(38) 



so the three-graviton vertex is consistent with the GR 
prediction at the 0.1% level. This bound is slightly worse, 
but comparable, to the one from lunar laser ranging, 
cq. (15). It should be stressed, however, that eq. (38) 
is really a test involving the radiative sector of GR, while 
eq. (15) only tests the conservative sector. For compari- 
son, observe that in the Standard Model the triple gauge 
boson couplings are measured to an accuracy of about 
3% [28]. 

For the double pulsar we find g(e) ~ 3.3. Since p, for 
the double pulsar is presently measured at the 1.4% level 
[8], we get a larger bound compared to eq. (38). However, 
further monitoring of this system is expected to bring the 
error on P, down to the 0.1% level. 



B. Binary coalescences at interferometers 

We now compare these results with what can be ex- 
pected from the detection of a binary coalescence at GW 
interferometers. In this case one can determine the phys- 
ical parameters of the inspiraling bodies, by performing 
matched filtering of theoretical waveform templates. In 
the matched filtering method any difference in the time 
behavior between the actual signal and the theoretical 
template model will eventually cause the two to go out 
of phase, with a consequent drop in the signal-to-noise 
ratio (SNR). The introduction of fa and fa affects the 
template, and in particular the accumulated phase 



2tt 



f(t)dt, 



(39) 



where f(t) is the time-varying frequency of the source, 
and the subscript min (max) denotes the values when 
the signal enters (leaves) the detector band-width. Thus 
in principle a detection of a GW signal from coalescing 
binaries could be translated into a measurement of the 
three- and four-graviton vertices. In this section we in- 
vestigate the accuracy of such a determination. 

With respect to the timing of binary pulsars, there are 
at least three important qualitative differences that af- 
fect the accuracy at which these systems can test the 



non-linearities of GR. First, coalescing compact bina- 
ries in the last stage of the coalescence reach values of 
v/c ~ 1/3, and are therefore much more relativistic than 
binary pulsars, which rather have v/c <~ 10~ 3 . Second, 
the leading Newtonian result for is of order (v/c) -5 
so it is much larger than one, and to get the phase with 
a precision A<fi <C 1, as needed by interferometers, all 
the corrections at least up to 0(v 6 /c 6 ) to the Newtonian 
result must be included, so higher-order corrections are 
important even if they are numerically small relative to 
the leading term. In other words, even if PN corrections 
are suppressed by powers of v/c with respect to the lead- 
ing term, they can be probed up to high order because 
what matters for GW interferometers is the overall value 
of the PN corrections to the phase, and not their value 
relative to the large Newtonian term. These two con- 
siderations should suggest that interferometers are much 
more sensitive than pulsar timing to the non-linearities 
of GR. 

On the other hand, for binary pulsars we can measure 
not only the decay of the orbital period due to GW emis- 
sion, but also several other Keplerian observables, that 
provide a determination of the geometry of the orbit, as 
well as post-Kcplerian observables, such as the periastron 
shift and the Einstein time delay, which fix the masses 
of the stars in the binary system. This is not the case 
for the detection of coalescences at interferometers. With 
interferometers the parameters that determine the wave- 
form, such as the masses and spins of the stars, must 
be determined from the phase of the GW itself, and one 
must then carefully investigate the degeneracies between 
the determination of fa (or of fa) and the determination 
of the masses and spins of the stars. This effect clearly 
goes in the direction of degrading the accuracy of param- 
eter reconstruction at GW interferometers, with respect 
to binary pulsar timing, so in the end it is not obvious 
a priori which of the two, GW interferometers or binary 
pulsar timing, is more sensitive to the non-linearities of 
GR. This question is answered in what follows. 

Repeating with fa 7^ the standard computation of 
the orbital phase for a circularized orbit, we find that 
fa modifies the orbital phase (f>(t) already at OPN (i.e. 
Newtonian) level, where to linear order in fa we get 



,0PN 



Q 5/8 



(1 + Kfa) 



with 60 = —5/2, and 6 is defined as 



9 = 



y{t c - 1) 
5GM 



(l - hfa) , 



(40) 



(41) 



where t c is the time of coalescence. Combining the factors 
v and M which enter in the definition of with the 
explicit factor l/v in cq. (40) we recover the well-know 
result that the Newtonian phase depends on the masses 
of the stars only through the chirp mass M c = v 3 / 5 M. 
From eq. (40) we immediately understand the crucial role 
that degeneracies have for interferometers. In fact, since 
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M c is determined from eq. (40) itself, using only the OPN 
phase (40) it is impossible to detect the deviation from 
the prediction of GR induced by /3 3 . A non-zero value of 
/3 3 would simply induce an error in the determination of 
M c . 

The same happens at 1PN level. In fact, at 1PN order 
and with (3% ^ the phase has the general form 



ilPN 



5/8 



where, as before, bo = 
portional to /3 3 , while 



(l + fe /3 3 ) + aiM(l+&i/3 3 )e- 1 ' /4 ] , 

(42) 

-5/2 is the OPN correction pro- 



, . 3715 55 
aiM = 8064 + 96" 



(43) 



is the 1PN GR prediction [17, 46], and b\ (which is pos- 
sibly ^-dependent) parametrizes the 1PN correction due 
to /3 3 . (For simplicity, we only wrote explicitly the term 
linear in /3 3 since |/3 3 | is much smaller than one, but all 
our considerations below can be trivially generalized to 
terms quadratic in (3 3 , just by allowing the function b\{u) 
to depend also on f3 3 ). In general we expect b\ to be O(l), 
and we will see below that for our purposes this estimate 
is sufficient. 

Using M c and v as independent mass variables, in place 
of mi and m 2 , we see that while the effect of (3 3 on the 
OPN phase can be reabsorbed into M c , its effect on the 
1PN phase can be reabsorbed into a rescaling of v. Ob- 
serve that, in the detection of a single coalescence event, 
GW interferometers do not measure the functional de- 
pendence of v of the 1PN phase, but only its numerical 
value for the actual value of v of that binary system, 
so we cannot infer the presence of a term proportional 
to from the fact that it changes the functional form 
of the ^-dependence from the one obtained by eq. (43). 
Thus, even at 1PN order, it is impossible to detect the 
deviations from GR induced by f3 3 . A non-zero (3 3 would 
simply induce an error on the determination of M c and 
u, i.e. on the masses of the two stars. 

We then examine the situation at 1.5PN order. Let us 
at first neglect the spin of the two stars. Then the 1.5PN 
phase with [3 3 ^ has the generic form 



<1.5PN 



5/S 



(l + 6 /3 3 ) + aiM(l + 6i/3 3 )e- 1/4 



+ a 2 (i + fo 2 /? 3 )e- 3 / 8 



(44) 



where a 2 = — 37r/4 is the 1.5PN GR prediction and 6 2 
is the (possibly {/-dependent) 1.5PN correction due to 
(3s. Again, we will not need its exact value, and we will 
simply make the natural assumption that it is O(l). 

However, for the purpose of determining /3 3 , neglecting 
spin is not correct. Indeed, at 1.5PN order the spin of 
the bodies enters through the spin-orbit coupling, and 
the evolution of the GW frequency / with time is given 



by [49] 
df _ 96 

di ~ y 



^8/3 M 5/3 / ll/3 1 _ 2 J_ ai[v)x + (4 ^ _ fa s)x m 

' (45)" 
where x — (irMf) 2 ' 3 , while /3ls describes the spin-orbit 
coupling and is given by 



As 



if 

12 ^ 



113 



M 2 



7bv 



(46) 



where L is the orbital angular momentum, \ a = ^o/ m o 
and S a is the spin of the a-th body. In principle /?ls 
evolves with time because of the precession of L, Si and 
S 2 . However, it turns out that in practice it is almost 
conserved, and can be treated as a constant [48]. Inte- 
gration of eq. (45) then shows that, in the 1.5PN phase, 
/3ls is exactly degenerate with /3 3 in eq. (44). Further- 
more, observe that /?ls, depending on the spin config- 
uration, can reach a maximum value of about 8.5 [48] 
(and its maximum value remains large even in the limit 
v — > 0), while /3 3 is already bound by laser ranging at 
the level of 2 x 10~ 4 and by pulsar timing at the level 
of 10~ 3 (which tests the radiative sector, as do GW in- 
terferometers). Thus, the effect of /3 3 at 1.5PN is simply 
reabsorbed into a (very small) shift of /3ls • 

At 2PN order (3 3 is degenerate with the parameter a 
that describes the spin-spin interaction 

a = ^[721(L-Xi)- 247 X i-X 2 (L-x 2 )] 



48 



96 



a=l 



2 
a 

M 2 



719(L- X J 2 - 233 X ^ • (47) 



The first term is the one which is usually quoted in the 
literature, first computed in [49] (see also [50, 51]). The 
term in the second line, computed recently in [52], is 
however of the same order, and must be included. 

The first term is proportional to v, and reaches a max- 
imum value cr max (^) ~ IOia In a coalescence with very 
small value of this term is therefore suppressed; e.g. 
in an extreme mass-ratio inspiral (EMRI) event at LISA 
where a BH of mass mi — IOMq falls into a supermas- 
sive BH with to 2 = 10 6 M Q , one has v = 10~ 5 and the 
term in the first line has a maximum value <~ 10~ 4 . If 
this standard term gave the full answer, a value of /3 3 
in excess of this value could therefore give an effect that 
cannot be ascribed to a. However, the presence of the 
new term recently computed in [52] spoils this reasoning, 
since it is not proportional to v. The conclusion is that, 
just as with /?sl at 1.5PN order, the effect of /3 3 at 2PN 
order is just reabsorbed into a small redefiniton of cr, and 
therefore simply induces an error in the reconstruction of 
the spin configuration (observe also that fixing /3ls does 
not allow us to fix the spin combinations that appear in 
a.) 

One could in principle investigate the effect of /3 3 on 
higher-order coefficients, such as the 2.5PN term -05 and 
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(1 + 04.) 



(a) 



(b) 



imation, in which only the harmonic at twice the source 
frequency is retained. Higher-order harmonics however 
break degeneracies between various parameters in the 
template [54], and it would be interesting to investigate 
whether their inclusion in the analysis allows one to put 
a bound on 03 and 04 from binary coalescences. 



FIG. 4: The diagrams that contribute to the conservative 
dynamics, which are affected by a modification of the four- 
graviton vertex. 

* ► 7 

{1 + l3i ]x3ULSULs 



FIG. 5: The diagram that contributes to the radiative dynam- 
ics, which is affected by a modification of the four-graviton 
vertex. 

the 3PN term ipQ in eq. (10), which at LISA can be mea- 
sured with a precision of order 1CT 2 [24]. Unfortunately, 
it is difficult to translate them into bounds on 03 because 
at 2.5PN order one finds that 03 is degenerate with a dif- 
ferent combination of spin and orbital variables, which is 
not fixed by the 1.5PN spin-orbit term (see Table 11 of 
ref. [53]), while at 3PN order the spin contribution is not 
yet known. The conclusion is therefore that interferome- 
ters cannot measure the three- and higher-order graviton 
vertex, since the effect of a modified vertex is simply re- 
absorbed into the determination of the masses and spin of 
the binary system. The conclusion that interferometers 
are not competitive with pulsar timing for measuring de- 
viations from GR was also reached in Ref. [10], although 
in a different context. In fact, Ref. [10] was concerned 
with multiscalar-tensor theories, whose leading-order ef- 
fect is the introduction of a term corresponding to dipolc 
radiation (a "minus one"-PN term). 

We now examine what can be said about the four- 
graviton vertex, parametrized by 04. In the conservative 
part of the Lagrangian 04 contributes through the dia- 
grams of Fig. 4. However, these contributions only affect 
the equations of motion at 2PN order, so there is no hope 
to see them in solar system experiments, where the ve- 
locities at play are very small. For the same reason, no 
significant bound can be obtained from binary pulsars; 
from a simple order of magnitude estimate we find that 
the Hulse- Taylor pulsar can only give a limit 04 < O(10), 
which is not significant. 

At GW interferometers, 04 enters into the phase for the 
first time at 1PN order, through the diagram in Fig. 5. 
However, it suffers exactly of the same degeneracy is- 
sues as 03, so it cannot be measured to any interesting 
accuracy at present or future interferometers, at least 
with the technique discussed here. Note however that in 
this paper we have worked in the restricted PN approx- 



V. CONCLUSION 

We have proposed to quantify the accuracy by which 
various experiments probe the non-linearities of GR, by 
translating their results into measurements of the non- 
abelian vertices of the theory, such as the three-graviton 
vertex and the four-graviton vertex. This is similar in 
spirit to tests of the Standard Model of particle physics, 
where the non-abelian vertices involving three and four 
gauge bosons have been measured at LEP and at the 
Tevatron. 

We have shown that, at a phenomenological level, this 
can be done in a consistent and gauge-invariant man- 
ner, by introducing parameters 03 and 04 that quantify 
the deviations from the GR prediction of the three- and 
four-graviton vertices, respectively. We have found that, 
in the conservative sector of the theory, i.e. as long as one 
neglects the emission of gravitational radiation at infin- 
ity, the introduction of 03 at 1PN order is phenomeno- 
logically equivalent to the introduction of a parameter 
/?ppn = 1 + 03 in the parametrized PN formalism. Strong 
bounds on 03 therefore come from solar system experi- 
ments, and most notably from lunar laser ranging, that 
provides a measurement at the 0.02% level. 

The modification of the three-graviton vertex however 
also affects the radiative sector of the theory, and we have 
found that the timing of the Hulse- Taylor pulsar gives a 
bound on 03 at the 0.1% level, not far from the one ob- 
tained from lunar laser ranging. Conceptually, however, 
the two bounds have different meanings, since lunar laser 
ranging only probes the conservative sector of the theory, 
while pulsar timing is also sensitive to the radiative sec- 
tor. 

We have then studied the results that could be ob- 
tained from the detection of coalescences at interferome- 
ters, and we have found that, even if 03 already modifies 
the GW phase at the Newtonian level and 04 at 1PN 
order, their effect can always be reabsorbed into other 
parameters in the template, such as the mass and spin 
of the two bodies so, rather than detecting a deviation 
from the GR prediction, one would simply make a small 
error in the estimation of these parameters. 
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At the quantum level, if a mass is generated, it will be 
power divergent with the cutoff, and will have to be fine 
tuned order by order in perturbation theory. However, 
the fact that quantum divergences have to be subtracted 
order by order is a generic problem of the standard quan- 
tum extension of GR, independently of fa. In any case, 
again, our approach is purely phenomenological, and the 
introduction of /3 3 is simply a tool for tracking a specific 
contribution to the computation. 



